Analysis of surface waves generated on subwavelength-structured silver films 
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Using transmission electron microscopy (TEM) to analyse the physical- chemical surface proper- 
ties of subwavlength structured silver films and finite- difference time-domain (FDTD) numerical 
simulations of the optical response of these structures to plane- wave excitation, we report on the 
origin and nature of the persistent surface waves generated by a single slit-groove motif and recently 
measured [H by far-field optical interferometry. The surface analysis shows that the silver films 
are free of detectable oxide or sulfide contaminants, and the numerical simulations show very good 
agreement with the results previously reported. 
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I. INTRODUCTION 

The optical response of subwavelength-structured metallic films has enjoyed a resurgence of interest in the past few 
years due to the quest for an all-optical solution to the inexorable drive for ever-smaller, denser integrated devices 
operating at ever-higher band width 2]. Two basic questions motivate research in this field: how to confine micron- 
sized light waves to subwavelength dimensions 3] and how to transmit this light without unacceptable loss over at 
least tens of microns 0, • Although surface waves and in particular periodic arrays of surface plasmon polaritons 
(SPPs) have received a great deal of attention as promising vehicles for subwavelength light confinement and transport, 
detailed understanding of their generation and early time evolution (within the first few optical cycles) in and on real 
metal films is still not complete 0, 0, H, ^] • Recent measurements ji[ of far-field interference fringes arising from 
surface wave generation in single slit-groove structures on silver films have characterized the amplitude, wavelength 
and phase of the surface waves. After a rapid amplitude decrease within the first three microns from the generating 
groove, waves persisting with near-constant amplitude over tens of microns were observed. Such long-range transport 
is the signature of a "guided mode" SPP, but the measured wavelength was found to be markedly shorter than the 
expectation from conventional theory One possible reason advanced for the disparity between experiment and 
theory was the presence of an oxide or sulfide dielectric layer on the silver surface, and it has been suggested 12] that 
an 11 nm layer of silver sulfide would bring experiment and theory into agreement. We report here the results of two 
investigations: one experimental, into the physical-chemical surface properties of the silver structures with the intent 
of detecting the presence of a dielectric layer and the other theoretical, into the calculated optical response of the 
untarnished silver slit-groove structures using the FDTD technique to numerically solve Maxwell's equations. These 
studies shown no evidence of oxide or sulfide layers on the silver surface and the numerical solutions to Maxwell's 
equations show good agreement with measurements reported in |1|] . 
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FIG. 1: Cross-sectional bright-field through- focal TEM micrograph of the Ag layer on a fused silica substrate. Labels refer to: 
fused silica substrate A, silver layer B, capping gold layer C, and capping platinum layer D. The darkened features in the Ag 
layer (B) are dislocations and grain boundaries. 

II. TEM ANALYSIS OF STRUCTURED SILVER SURFACES 

Three typical structured samples were chosen for examination with fabrication dates of about 12 months, 6 months, 
and 1 week from the date of the TEM analysis. The structures dating from 6 months and 12 months were part of the 
series of structures actually used in the previous reports 0, . The subwavelength structures consist of a 400 nm 
silver film evaporated onto a fused silica substrate (Corning 7980 UV grade® 25 mm square, 1 mm thick optically 
polished on both sides to a roughness of no more than 0.7 nm). The structured substrates are stored in Fluoroware® 
sample holders, 25 mm diameter and 1 mm in depth at the center. 

A. Methods 

Electron transparent sections for cross-sectional transmission electron microscopy (TEM) examination were pre- 
pared by sample thinning to electron transparency using standard focused Ga+-ion beam (FIB) milling procedures 
[l^ in a FEI 200 FIB workstation and placed on a holey carbon support. The TEM characterization was performed 
using a Philips CM300 Schottky field emission gun (FEGTEM) microscope operating at 300 kV. The field emission 
electron source is ideal for applications requiring high coherency, high brightness at high magnification, or small fo- 
cused probes of selected areas. The FEGTEM has a point resolution of 0.2 nm and an information limit of 0.12 nm. 
The minimum focused electron beam probe size is 0.3 nm. In dark field images, one or more diffracted beams are 
allowed to pass the objective aperture while the direct straight-through beam is blocked by the aperture. In contrast 
to the direct beam, the diffracted beam interacts strongly with the specimen and selection of a particular diffracted 
beam allows better visual phase differentiation. 

B. Examination of silver (Ag) layers 

Typical regions of Ag layer that were thinned to electron transparency are shown at relatively low magnification in 
the bright field micrograph in Fig.^ An immediate distinction can be made between the fused silica substrate (A) 
and the Ag deposit (B). Protective capping layers of Au and Pt, applied at the time of analysis, are marked at (C) 
and (D) respectively. The latter are used to prevent any 'top down' ion damage of the cross-section during the ion 
beam thinning preparation. 

The white colored line marked by an arrow in Fig.Qand shown in more detail in Fig.Ois a band that contains very 
fine particles of Au that are formed during the initial stages of the deposition of this metal. 





FIG. 2: Cross-sectional bright-field through- focal 
TEM micrograph of the Ag layer (B) on the fused 
silica substrate (A). The protective capping layers 
of Au and Pt are marked at (C) and (D) respec- 
tively. 
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FIG. 3: A higher magnification of the upper sur- 
face showing the Au particle-containing protec- 
tion layer. The undulations and voids indicate 
a surface nonuniformity somewhat greater than 
that of the samples in general (± ^ 10 nm from 
400 nm nominal thickness) but the site was chosen 
to more vividly show the initial line of gold capping 
nanoparticles. 





FIG. 4: Cross-sectional bright-field through-focal TEM micrographs of the Ag layer on the fused silica substrate. The images 
were acquired at over- focus (left), in- focus (middle) and under- focus (right). 



The fine particles do not exhibit the same degree of absorption contrast as the bulk of the Au formed above them 
because the particle diameters are less than the thickness of the cross-sectional slice and thus do not extend through 
the full thickness of the TEM sample, as do the relatively coarse grains constituting the bulk of the Au protection 
layer. A higher magnification TEM micrograph of this region is shown for clarity in Fig.IHl It is clearly observed 
that the Au nanoparticles form a separate layer above the Ag deposit and any local undulations in the Ag layer are 
observed to be devoid of any oxide or sulfide layer. 

Fresnel contrast methods were used to examine the upper surface regions of the silver in more detail and a typical 
area is shown in Fig.^ Imaging of layers in cross-section, in which the electrons travel parallel to the interfaces 
between the two materials, always results in Fresnel fringes. Because the amplitude and phase changes that occur 
when an electron is scattered elastically are characteristic of the atomic number, there will be changes in the elastic 
scattering directly related to the form of the projected scattering potential when there is a composition change at an 
interface viewed in projection 15], analogous to a phase grating filtering of the scattered electronic wave function. 
Interface Fresnel effects can thus provide a signature on the form and magnitude of any compositional discontinuity 
that is present. The visibility of these fringes depends on the thickness of the specimen and on the defocus value of 
the microscope. The contrast of the fringes is higher with increasing defocus, at both positive and negative values 




FIG. 5: Cross-sectional bright-field through- focal TEM micrographs of the tarnished Ag layer on the fused silica substrate after 
more than two years exposed to ambient laboratory conditions. Insert: electron diffraction pattern of the tarnished Ag layer 
showing amorphous (noncrystalline) characteristics by means of broad, diffuse diffraction rings. 

from the minimum contrast condition. From Fig.^ the critical observation is in the absence of Fresnel effects at the 
surface of the metal as a function of focus conditions other than at the localized regions associated with the gold 
protective coating. The lack of Fresnel fringes at the upper surface indicates a lack of variation in the scattering 
potentials and hence of the chemistry or the composition at the interface. The lack of Fresnel signature demonstrates 
that silver surface is not covered by a detectable layer of sulfide or oxide since their presence would exhibit differences 
in scattering potential by comparison with the metal itself The lower detection limit for such compositional 
differences is a few tenths nanometer layer thickness. 

C. Discussion of the surface analysis 

The high resolution TEM analysis failed to detect the presence of oxide or sulfide coating on any of the structured 
silver samples. This result is in fact not surprising and is consistent with a simple estimate of the surface chemical 
impurity likely to be present. Since the structures are kept in air-tight sample holders, the total, static volume of air 
to which they are exposed during storage is 2000 mm^. The trace fractional concentrations of sulfur dioxide and 
hydrogen sulfide in ordinary laboratory air is less than 0.2 x 10~^ and therefore about 10^^ impurity molecules are 
available to react with the 625 mm^ silver surface. The average linear distance between impurity molecules, assuming 
all of them reacted with the surface, is therefore about 250 nm - a surface density below the detection limit of the 
FEGTEM instrument. 

For comparison, TEM was also conducted on a similar Ag structure left fully exposed to ambient laboratory 
environment for more than two years. The resulting micrograph and electron diffraction pattern are shown in Fig.ISl 
The tarnished Ag layer is observed to be uniform and the corresponding electron diffraction data shows a fully 
amorphous (noncrystalline) character evidenced by the broad diffuse rings. The detection of a relatively high intensity 
signal from a nominally nondiffracting material yields valuable information on not only the crystallinity of the material, 
but also the chemical composition. Comparing the electron diffraction data for the tarnished samples to that for the 
silver structures described earlier, we find no evidence of a similar noncrystalline material on the surface of the 
untarnished Ag shown in Fig.^ 

III. CALCULATION OF THE FAR-FIELD TRANSMISSION 

The optical response of structured metal surfaces is simulated using a finite-difference-time-domain (FDTD) ap- 
proach pTl. The subwavelength slit-groove structures l„iliij can be modelled in two dimensions, and we consider here 
the case of transverse magnetic TM polarization (with the i7-field transverse to the plane of incidence and parallel to 
the long axis of the slit-groove structures). The set-up for the simulations is shown in Fig.El The Maxwell equations 



far field detection screen 
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FIG. 6: Schematic structure and detection contour for FDTD simulations. The sht and the groove widths are 100 nm, the 
groove depth is 100 nm and the silver film thickness is 400 nm. Interference of waves directly propagating through the slit with 
waves originating from the groove and propagating to the slit as surface waves give rise to the simulated interference pattern 
shown in Fig.|3 



are thus, 
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where Ex^ Ey^ and Hz are the Cartesian components of the electric and magnetic fields, Jx and Jy are the Cartesian 
components of the current density J, e is the electric permittivity and jhq is the magnetic permeability of free space. 
In metallic regions of space £ is a complex valued, frequency-dependent function. Within the standard Drude model 
[l8j it is given as. 



where Sq is the electric permittivity of free space, Soo = £ {co ^ oo) is the dimensionless infinite frequency limit of the 
dielectric constant, ujp is the bulk plasmon frequency, and F is the damping rate. In order to properly describe the 
material dispersion of the metal structures we numerically fit the real and imaginary parts of the Drude dielectric 
constant in the form of Eq.O to the experimental data collected in 0| . The Drude parameters obtained in the 
wavelength regime ranging from 750 nm to 900 nm for silver are Soo = 3.2938, cOp = 1.3552 x 10^^ rad s~^, and 
F = 1.9944 X 10^^ rad s-\ which correspond to Re [e] = -33.9767 and Im [e] = 3.3621 at Aq = 852 nm. These fitted 
parameters are close to those determined in Q by ellipsometry, Re[e] = —33.27 and lm[e] = 1.31. Numerical results 
were not sensitive to this range of parameter variability. 

The frequency dependence of the dielectric constant Eq.[2l results in an additional equation to the Maxwell equa- 
tions (P), describing the time evolution of the current density J in the metallic region of space [20^ . 



where a = — F, (3 = eoco?, and e = £0^00 • In the surrounding free space e = Sq^ and a = = and hence J vanishes 

m- 

As a reference for the optical response of the Drude metal we also implement perfect electric conductor (PEC) 
boundary conditions, where the metal dielectric constant is set to negative infinity, and all electromagnetic field 
components are therefore strictly zero in metal regions. 

The light source is a plane wave, incident perpendicular to the metal surface and depending on time as, 

E,^,{t)=Eof{t)cosut, (4) 

where £^0 is the peak amplitude of the pulse, f {t) = sin^ (7rt/r) is the pulse envelope, and r is the pulse duration 
(/(t>r)=0. 

Time propagation of Eqs.dOlis performed by a leapfrogging technique [13 • In order to prevent nonphysical refiection 
of outgoing waves from the grid boundaries, we employ perfectly matched layer (PML) absorbing boundaries, within 
which Hz is split into two additive subcomponents [23 • To avoid accumulation of spurious electric charges at the 
ends of the excitation line and generate a pure plane wave with a well-defined incident wavelength, we embed the 
ends in the PML regions. We have tested this approach by direct comparison with the total field/scattered field 
technique 0. In all simulations the spatial grid ranges from x = —16.2 /im to 16.2 /im and y = —1.79 /im to 1.79 /im. 
Convergence is achieved with a spatial step size of Sx = Sy ^ 4nm and a temporal step size of St = Sx/ (2c), where c 
denotes the speed of light in vacuum. 

All simulations have been performed on distributed memory parallel computers at the National Energy Research 
Scientific Computing Center and San Diego Supercomputer Center. The parallel technique used in our simulations is 
described in detail in ^23.] . 



A. Numerical results 



Calculations of the intensity, / ~ (£^^ -\- Ey -\- H^), are performed along a box contour as shown in Fig. including 
a line parallel to the metal surface and located at the distance Ld- The collected data is averaged over time and the 
spatial coordinates for a range of slit-groove distances. Finally, the space- and time-averaged intensity is normalized 
to unit maximum. Our results converge with Ld > 1 jam and are invariant to the pulse duration (hence converged to 
the CW limit) with incident pulse durations r > 200 fs. A direct comparison of the experimental data with both the 
PEC and the Drude models is shown in Fig.[3 

Clearly, the Drude model agrees very well with the data, whereas the PEC model predicts oscillations of the 
intensity with a noticeably larger wavelength and smaller amplitude. In order to properly determine the wavelength 
of oscillations we use a cos function with an exponentially decreasing amplitude component and constant offset to fit 
the data shown in Fig.[3and extend both the FDTD simulation and the fitting function to a slit-groove distance of 
16 /im, well beyond the range of experimental data. We then take the Fourier transform of the fitted function over 
this extended range to obtain the associated power spectrum expressed as a wavelength distribution. The results are 
shown in Fig.IHl 

ht{x) = [Ai + A2 exp(A3x)] cosiA^x + A5) + As (5) 

The FDTD simulation and analytic fit for larger slit-groove distances are shown in Fig.|Hlfor the case of the Drude 
model. The fitting function Eq.l^l tracks reasonably well the FDTD results over the entire slit-groove distance range. 
Figureini shows the normalized power spectrum corresponding to Eq.ISl as a function of the wavelength. The two 
power spectrum plots of Fig. [HI provide the effective wavelength, Aeff, at which surface waves propagate as well as 
the distribution of modes around the peak. Note that the Drude model for silver shows a marked blue shift in peak 
wavelength and a noticeable broadening of the distribution compared to the perfect metal PEC model. The Drude 
model results in Aeff = 837.482 nm , whereas PEC model gives Aeff = 852.066 nm, very close to the free-space reference 
wavelength Aq = 852 nm. The effective surface index of refraction 

^eff T , (6) 

Aeff 

where Ainc is the incident wavelength (in the experiment Ainc = 852 nm), leads to the following values: neff (Drude) = 
1.0173 and neff (PEC) = 0.9999. 
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FIG. 7: Comparison of FDTD simulation and experiment. Green points are experimental data taken froml^. Blue curve is 
the FDTD result for the Drude model. Red curve is the FDTD result for the PEC model. 
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FIG. 8: Blue curve plots the same Drude model FDTD calcu- 
lation as in Fig.UJbut extended to 16 fim slit-groove distance. 
Red dashed curve plots the the analytic fitting function of 
Eq.El 
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FIG. 9: Fourier transform spectra of Eq.^for curves fitted 
to the Drude (blue) and PEC (red) results of FDTD simu- 
lations. Fitting parameters used in Eq.^for Drude (PEC) 
models are, Ai:0.20686 (0.061126), ^2:0.19597 (0.22834), 
^31-0.0005213 (-0.00046489), ^4:0.007374 (0.0075026), 
^5:3.3098 (2.0663), ^6:0.54336 (0.71776). 



IV. SUMMARY AND CONCLUSIONS 



Surface analysis by transmission electron microscopy of subwavelength structured silver films used to investigate 
their optical response 1, 10] showed no detectable evidence of material on the surface other than silver. The suggestion 
that an 11 nm sulfide layer may be present so as to bring the interference pattern calculated by the authors of Ref. j3| 
into agreement with experiment is therefore not confirmed. 

In contrast to the calculations reported by Ref. 0, numerical solution of Maxwell's equations reported here using 
an FDTD approach does show very good agreement with the fringe amplitude and wavelength over the slit-groove 
range of the data reported in ^ . Extrapolation of the FDTD simulations beyond the range of the measurements shows 
that the initially decreasing amplitude of the fringe settles to an oscillation with near-constant amplitude and fringe 



contrast. These features resemble the calculations of but the results reported here show much better agreement 
with experiment without the need to invoke an 11 nm silver sulfide layer. Fourier analysis of the FDTD simulations 
reveal that the most probable wavelength in the Fourier distribution is 837 nm, within a about 2 nanometers of the 
expected long-range SPP wavelength of 839 nm. The reasons for the discrepancies between the results of Ref. [T^ . 
which employ a different approach using rigorous-coupled-wave-analysis (RCWA) and those reported here have not 
yet been identified. 
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